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Abstract 

We propose a generalization of the random matrix theory following the basic prescription of 
the recently suggested concept of superstatistics. Spectral characteristics of systems with mixed 
regular-chaotic dynamics are expressed as weighted averages of the corresponding quantities in the 
standard theory assuming that the mean level spacing itself is a stochastic variable. We illustrate 
the method by calculating the level density, the nearest-neighbor-spacing distributions and the 
two- level correlation functions for system in transition from order to chaos. The calculated spacing 
distribution fits the resonance statistics of random binary networks obtained in a recent numerical 
experiment. 
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I. INTRODUCTION 



Random matrix theory (RMT) provides a suitable framework to describe quantal systems 
whose classical counterpart has a chaotic dynamics 

HQ. 

It models a chaotic system by 
an ensemble of random Hamiltonian matrices H that belong to one of the three universal 
classes, namely the Gaussian orthogonal, unitary and symplectic ensembles (GOE, GUE 
and GSE). The theory is based on two main assumptions: (i) the matrix elements are 
independent identically-distributed random variables, and (ii) their distribution is invariant 
under unitary transformations. These lead to a Gaussian probability density distribution for 
the matrix elements, P (H) oc exp [— 77 Tr (H^H)]. With these assumptions, RMT presents 
a satisfactory description for numerous chaotic systems. On the other hand, there are 
elaborate theoretical arguments by Berry and Tabor which are supported by several 
numerical calculations, that the nearest-neighbor-spacing (NNS) distribution of classically 
integrable systems should have a Poisson distribution exp(— s), although exceptions exist. 

For most systems, however, the phase space is partitioned into regular and chaotic do- 
mains. These systems are known as mixed systems. Attempts to generalize RMT to describe 
such mixed systems are numerous; for a review please see 0. Most of these attempts are 
based on constructing ensembles of random matrices whose elements are independent but 
not identically distributed. Thus, the resulting expressions are not invariant under base 
transformation. The first work in this direction is due to Rosenzweig and Porter [5]. They 
model the Hamiltonian of the mixed system by a superposition of a diagonal matrix of 
random elements having the same variance and a matrix drawn from a GOE. Therefore, 
the variances of the diagonal elements of the total Hamiltonian are different from those of 
the off-diagonal ones, unlike the GOE Hamiltonian in which the variances of diagonal el- 
ements are twice those of the off-diagonal ones. Hussein and Pato |6[ used the maximum 
entropy principle to construct "deformed" random- matrix ensembles by imposing different 
constraints for the diagonal and off-diagonal elements. This approach has been successfully 
applied to the case of metal- insulator transition 71. A recent review of the deformed en- 
semble is given in [8j. Ensembles of band random matrices, whose entries are equal to zero 
outside a band of limited width along the principal diagonal, have often been used to model 
mixed systems ODD However, so far in the literature, there is no rigorous statistical 
description for the transition from integrability to chaos. The field remains open for new 
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proposals. 

The past decade has witnessed a considerable interest devoted to the possible general- 
ization of statistical mechanics. Much work in this direction followed Tsallis seminal paper 
~l 

Tsallis introduced a non-extensive entropy, which depends on a positive parameter q 
known as the entropic index. The standard Shannon entropy is recovered for q = 1. Ap- 
plications of the Tsallis formalism covered a wide class of phenomena; for a review please 



see, e.g. 



12| . Recently, the formalism has been ar 



led to include systems with mixed 



regular-chaotic dynamics in the framework of RMT [13J, [14j, [15|, [16J, [17|, [18( . This is done by 
extremizing Tsallis' non-extensive entropy, rather than Shannon's, but again subject to the 
same constraints of normalization and existence of the expectation value of Tr(if^). The 
latter constraint preserves base invariance. The first attempt in this direction is probably 
due to Evans and Michael Q]. Toscano et al. 14 1 constructed non-Gaussian ensemble 
by minimizing Tsallis' entropy and obtained expressions for the level densities and spacing 
distributions for mixed systems belonging to the orthogonal-symmetry universality class. 
Bertuola et al. expressed the spectral fluctuation in the subextensive regime in terms 
of the gap function, which measures the probability of an eigenvalue-free segment in the 
spectrum. A slightly different application of non-extensive statistical mechanics to RMT is 
due to Nobre et al. 16]. The nearest-neighbor-spacing (NNS) distributions obtained in this 
approach decays as a power-law for large spacings. Such anomalous distributions can hardly 
be used to interpolate between nearly-regular systems which have almost exponential NNS 
distributions and nearly-chaotic ones whose distributions behave at large spacing as Gaus- 
sians. Moreover, the constraints of normalization and existence of an expectation value for 
Ti(WH) set up an upper limit for the entropic index q beyond which the involved integrals 
diverge. This restricts the validity of the non-extensive RMT to a limited range near the 
chaotic phase jl?! I3| • 

Another extension of statistical mechanics is provided by the formalism of superstatistics 
(statistics of a statistics), recently proposed by Beck and Cohen (lflj . Superstatistics arises 
as weighted averages of ordinary statistics (the Boltzmann factor) due to fluctuations of 
one or more intensive parameter (e.g. the inverse temperature). It includes Tsallis' non- 
extensive statistics, for q > 1, as a special case in which the inverse temperature has a 
^-distributions. With other distributions of the intensive parameters, one comes to other 
more general superstatistics. Generalized entropies, which are analogous to the Tsallis 
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entropy, can be defined for these general superstatistics 
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This formalism has 



)een elaborated and applied successfully to a wide variety of physical problems, e.g., in 



2a |2J, |2jt |2a I23, |28l |2 



In a previous paper |3l|, the concept of superstatistics was applied to model a mixed 
system within the framework of RMT. The joint matrix element distribution was represented 
as an average over exp [— r/Tr (iff if] ] with respect to the parameter r\. An expression for 
the eigenvalue distributions was deduced. Explicit analytical results were obtained for the 
special case of two-dimensional random matrix ensembles. Different choices of parameter 
distribution, which had been studied in Beck and Cohen's paper [3] were considered. These 
distributions essentially led to equivalent results for the level density and NNS distributions. 
The present paper is essentially an extension of the superstatistical approach of Ref. 31 1 
to random-matrix ensembles of arbitrary dimension. The distribution of local mean level 
densities is estimated by applying the principle of maximum entropy, as done by Sattin 
[22]]. In Section 2 we briefly review the superstatistics concept and introduce the necessary 
generalization required to express the characteristics of the spectrum of a mixed system 
into an ensemble of chaotic spectra with different local mean level density. The evolution 
of the eigenvalue distribution during the stochastic transition induced by increasing the 
local-density fluctuations is considered in Section 3. The corresponding NNS distributions 
are obtained in Section 4 for systems in which the time-reversal symmetry is conserved or 
violated. Section 5 considers the two-level correlation functions. The conclusion of this work 
is formulated in Section 6. 



II. FORMALISM 

A. Superstatistics and RMT 

To start with, we briefly review the superstatistics concept as introduced by Beck and 
Cohen Consider a non-equilibrium system with spatiotemporal fluctuations of the 

inverse temperature (3. Locally, i.e. in spatial regions (cells) where j3 is approximately 
constant, the system may be described by a canonical ensemble in which the distribution 
function is given by the Boltzmann factor e~ f3E , where E is an effective energy in each 
cell. In the long-term run, the system is described by an average over the fluctuating f3. 
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The system is thus characterized by a convolution of two statistics, and hence the name 
"superstatistics" . One statistics is given by the Boltzmann factor and the other one by the 
probability distribution f{(3) of f3 in the various cells. One obtains Tsallis' statistics when f3 
has a x 2 distribution, but this is not the only possible choice. Beck and Cohe n g ive several 



possible examples of functions which are possible candidates for f(f3). Sattin [27jj suggested 
that, lacking any further information, the most probable realization of f{fi) will be the one 
that maximizes the Shannon entropy. Namely this version of superstatistics formalism will 
now be applied to RMT. 

Gaussian random-matrix ensembles have several common features with the canonical 
ensembles. In RMT, the square of a matrix element plays the role of energy of a molecule 
in a gas. When the matrix elements are statistically identical, one expects them to become 
distributed as the Boltzmann's. One obtains a Gaussian probability density distribution of 
the matrix elements 

P (H) oc exp [-rjTi (H^H)] (1) 

by extremizing the Shannon entropy ]} 32 ] subjected to the constraints of normalization 
and existence of the expectation value of Tr(if^if). The quantity Tr (if^il) plays the role 
of the effective energy of the system, while the role of the inverse temperature (3 is played 
by 7], being twice the inverse of the matrix-element variance. 

Our main assumption is that Beck and Cohen's superstatistics provides a suitable de- 
scription for systems with mixed regular-chaotic dynamics. We consider the spectrum of 
a mixed system as made up of many smaller cells that are temporarily in a chaotic phase. 
Each cell is large enough to obey the statistical requirements of RMT but has a different 
distribution parameter rj associated with it, according to a probability density f(rj). Con- 
sequently, the super statistical random-matrix ensemble that describes the mixed system is 
a mixture of Gaussian ensembles. Its matrix-element joint probability density distributions 
obtained by integrating distributions of the form in Eq. (1) over all positive values of rj with 
a statistical weight /(r/), 

where Z{rj) = J exp [— r/Tr (iffiif)] drj. Here we use the "B-type superstatistics" The 
distribution in Eq. (2) is isotropic in the matrix-element space. Relations analogous to 
Eq. (1) can also be written for the joint distribution of eigenvalues as well as any other 



statistic that is obtained from it by integration over some of the eigenvalues, such as the 
nearest-neighbor-spacing distribution and the level number variance. The distribution f(rj) 
has to be normalizable, to have at least a finite first moment 

(v) = / f(v)vdv, (3) 
Jo 

and to be reduces a delta function as the system becomes fully chaotic. 

The random-matrix distribution in Eq. (2) is invariant under base transformation be- 
cause it depends on the Hamiltonian matrix elements through the base-invariant quantity 
Tr (H^H). Factorization into products of individual element distributions is lost here, un- 
like in the distribution functions of the standard RMT and most of its generalizations for 
mixed systems. The matrix elements are no more statistically independent. This handi- 
caps one in carrying numerical calculations by the random-number generation of ensembles 
and forces one to resort to artificial methods as done in Q|. Base invariance makes the 
proposed random-matrix formalism unsuitable for description of nearly integrable systems. 
These systems are often described by an ensemble of diagonal matrices in a presumably 
fixed basis. For this reason we expect the present superstatistical approach to describe only 
the final stages of the stochastic transition. The base invariant theory in the proposed form 
does not address the important problem of symmetry breaking in a chaotic system, where 
the initial state is modelled by a block diagonal matrix with m blocks, each of which is 
a GOE m. This problem is well described using deformed random- matrix ensembles as 

n 

in p or phenomenologically by considering the corresponding spectra as superpositions of 
independent sub-spectra, each represented by a GOE |33j . 

The physics behind the proposed superstatistical generalization of RMT is the following. 
The eigenstates of a chaotic system are extended and cover the whole domain of classically 
permitted motion randomly, but uniformly. They overlap substantially, as manifested by 
level repulsion. There are no preferred eigenstate; the states are statistically equivalent. As a 
result, the matrix elements of the Hamiltonian in any basis are independently but identically 
distributed, which leads to the Wigner-Dyson statistics. Coming out of the chaotic phase, 
the extended eigenstates become less and less homogeneous in space. Different eigenstates 
become localized in different places and the matrix elements that couple different pairs are 
no more statistically equal. The matrix elements will no more have the same variance; one 
has to allow each of them to have its own variance. But this will dramatically increase 
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the number of parameters of the theory. The proposed superstatistical approach solves this 
problem by treating all of the matrix elements as having a common variance, not fixed but 
fluctuating. 

B. Eigenvalue distribution 

The matrix-element distribution is not directly useful in obtaining numerical results con- 
cerning energy-level statistics such as the nearest-neighbor spacing distribution, the two- 
point correlation function, the spectral rigidity, and the level-number variance. These quan- 
tities are presumably obtainable from the eigenvalue distribution. From (1), it is a simple 
matter to set up the eigenvalue distribution of a Gaussian ensemble. With H = U~ X EU, 
where U is the global unitary group, we introduce the elements of the diagonal matrix of 
eigenvalues E = diag(£i, • • ■ , E^) of the eigenvalues and the independent elements of U as 
new variables. Then the volume element (4) has the form 



where A N (E) = Y\ n>m {E n — E m ) is the Vandermonde determinant and dfi(U) the invariant 
Haar measure of the unitary group Q, I^J. Here (3 — 1,2 and 4 for GOE, GUE and GSE, 
respectively. The probability density Pp{H) is invariant under arbitrary rotations in the 
matrix space. Integrating over U yields the joint probability density of eigenvalues in the 
form 



where Pp '(rj, Ei, ■ ■ ■ , E N ) is the eigenvalue distribution of the corresponding Gaussian en- 
semble, which is given by 



where Cp is a normalization constant. Similar relations can be obtained for any statistic 
ap(Ei, ■ ■ • , Ef.), with k < N, that can be obtained from Pp(E 1: ■ • • , E N ) by integration over 
the eigenvalues E k+ i, • • • , E N . 

In practice, one has a spectrum consisting of a series of levels {Ei} , and is interested in 
their fluctuation properties. In order to bypass the effect of the level density variation, one 



dH = \ A N (E)f dEdfi(U), 



(4) 
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introduces the so called "unfolded spectrum" where Si = Ei/D and D is the local mean 
level spacing. Thus, the mean level density of the unfolded spectrum is unity. On the other 
hand, the energy scale for a Gaussian random-matrix ensemble is defined by the parameter 
7]. The mean level spacing may be expressed as 

D=^=, (7) 

where c is a constant depending on the size of the ensemble. Therefore, although the 
parameter rj is the basic parameter of RMT, it is more convenient for practical purposed to 
consider the local mean spacing D itself instead of rj as the fluctuating variable for which 
superstatistics has to be established. 

The new framework of RMT provided by superstatistics should now be clear. The local 
mean spacing D is no longer a fixed parameter but it is a stochastic variable with probability 
distribution f(D). Instead, the the observed mean level spacing is just its expectation 
value. The fluctuation of the local mean spacing is due to the correlation of the matrix 
elements which disappears for chaotic systems. In the absence of these fluctuations, f(D) = 
5(D — 1) and we obtain the standard RMT. Within the superstatistics framework, we can 
express any statistic a(E) of a mixed system that can in principle be obtained from the 
joint eigenvalue distribution by integration over some of the eigenvalues, in terms of the 
corresponding statistic a^ G \E,D) for a Gaussian random ensemble. The superstatistical 
generalization is given by 



The remaining task of superstatistics is the computation of the distribution f(D). 
C. Evaluation of the local- mean-spacing distribution 

Following Sattin j^J, we use the principle of maximum entropy to evaluate the distri- 
bution f(D). Lacking a detailed information about the mechanism causing the deviation 
from the prediction of RMT, the most probable realization of f(D) will be the one that 
extremizes the Shannon entropy 




(8) 




(9) 
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with the following constraints: 

Constraint 1. The major parameter of RMT is rj defined in Eq. (1). Superstatistics 
was introduced in Eq. (2) by allowing r\ to fluctuate around a fixed mean value (77). This 
implies, in the light of Eq. (7), the existence of the mean inverse square of D, 

(D- 2 ) = J o f(D)—dD. (10) 

Constraint 2. The fluctuation properties are usually defined for unfolded spectra, which 
have a unit mean level spacing. We thus require 



f(D)DdD = 1. (11) 
Therefore, the most probable f(D) extremizes the functional 

/oo poo poo 1 

f(D)lnf(D)dD-\ 1 j f(D)DdD-\ 2 J^ f(D)—dD (12) 



where Ai and A2 are Lagrange multipliers. As a result, we obtain 



f(D) = Cexp 



—a 



[2D DgM 



(13) 



where a and D are parameters, which can be expressed in terms of the Lagrange multipliers 
Ai and A2, and C is a normalization constant. We determine D and C by using Eqs. (10) 
and (11) as 

G30(a 3 |0,i,l) , x 

D —^ti' (14) 



and 



D Q GH (o»| 0,1,1)" 1 " J 



Here (x \ b\, 6 2 , b 2 ) is a Meijer's G-function defined in the Appendix. 



III. LEVEL DENSITY 



The density of states can be obtained from the joint eigenvalue distribution directly by 
integration 

p{E) =nJ -J P P {E, Ei,---, E N )dE 2 ■ ■ ■ dE N . (16) 
For a Gaussian ensemble, simple arguments Q, Q lead to Wigner's semi-circle law 



. ^JRl-E 2 , for \E\ < R 
Pge(E,D)={ ^ V , (17) 

0, for \E\ > R 



where D is the mean level spacing, while the prefactor is chosen so that p GE (E) satisfies the 
normalization condition 

/oo 
p GE (E)dE = N. (18) 
-oo 

We determine the parameter R Q by requiring that the mean level density is 1/D so that 

1 Z" 00 1 

N J J PGE{E)]2 dE = D- (19) 

This condition yields 

Ro = ^D. (20) 

Substituting (17) into (8) we obtain the following expression for the level density of the 
superstatistical ensemble 

/■37r 2 |E|/(16JV) 

p SE (E,a) = f{D,a)p GE (E,D)dD. (21) 

Jo 

We could not solve this integral analytically. We evaluated it numerically for different values 
of a. The results of calculation are shown in Fig. 1. The figure shows that the level density 
is symmetric with respect to E — for all values of a and has a pronounced peak at the 
origin. However, the behavior of the level density for finite a is quite distinct from the 
semicircular law. It has a long tail whose shape and decay rate both depend on the choice 
the parameter distribution f(D). This behavior is similar to that of the level density of 
mixed system modelled by a deformed random- matrix ensemble 34]. 



IV. NEAREST-NEIGHBOR-SPACING DISTRIBUTION 

The NNS distribution is probably the most popular characteristic used in the analysis 
of level statistics. In principle, it can be calculated once the joint-eigenvalue distribution is 
known. The super statistics generalization of NNS distribution for an ensemble belonging to 
a given symmetry class is obtained by substituting the NNS distribution of the corresponding 
Gaussian ensemble P GE (s, D) for a^ G \E, D) in (7) and integrating over the local mean level 
spacing D 

POO 

P SE (s) = / f{D)P GE (s,D)dD. (22) 
Jo 

Till now, no analytical expression for the NNS distribution could be derived from RMT. 
What we know is that this distribution is very well approximated by the Wigner surmise 
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l|. We shall obtain superstatistics for NNS distribution for systems with orthogonal and 
unitary symmetries by assuming that the corresponding Gaussian ensembles have Wigner 
distributions for the nearest-neighbor spacings. 

Equation (22) yields the following relation between the second moment (D 2 ) of the local- 
spacing distribution f(D) and the second moment (s 2 ) of the spacing distribution Pse(s): 

= (23) 

where (s 2 ) GE is the mean square spacing for the corresponding Gaussian ensemble. For the 
distribution in Eq. (13), one obtains 

= GgO (a 3 |0; l )l)G 30 (a 3 |0; 3 )2) 

^ ' K(a3| ,l,|)] 2 

Using the asymptotic behavior of the G-function, we find that (D 2 ) — > 1 as a — > oo, while 
(D 2 ) = 2 (as for the Poisson distribution) when a = 0. For practical purposes, the expression 
in Eq.(24) can be approximated with a sufficient accuracy by (D 2 ) « 1 + 1/(1 + 4.121a). 
Thus, given an experimental or numerical-experimental NNS distibution, one can evaluate 
the quantity (s 2 ) and estimate the corresponding value of the parameter a by means of the 
following approximate relation 

"""■^(PRkr (25) 

A. Orthogonal ensembles 

Systems with spin-rotation and time-reversal invariance belong to the orthogonal sym- 
metry class of RMT. Chaotic systems of this class are modeled by GOE for which NNS is 
well approximated by the Wigner surmise 

W^) = ^exp(-^ 2 ). (26) 

We now apply superstatistics to derive the corresponding NNS distribution assuming that 
the local mean spacing distribution f(D) is given by Eq. (13). Substituting (26) into (22), 
we obtain 



not 2 ( 1XOL 2 

Psoe(s, a) = — — — ___sGg2 a 3 + 



2^G30( a 3| ,I,l)^ 03 V "4D§" 
11 



2 



,0,0 , (27) 



where D is given by (14), while the suffix SOE stand for Superstatistical Orthogonal En- 
semble. 

Because of the difficulties of calculating Gq'® (z b 2 , b 3 ) at large values of z, we use (say 
for z > 100) the large z asymptotic formula given in the Appendix to obtain 



exp 

7T 

PsoE{s,a) « -s 



-3a ( A 7i + fr-i 



2 /-I 1 7T5 2 

1 ^ 4a 



(28) 



which clearly tends to the Wigner surmise for the GOE as a approaches infinity. This 
formula provides a reasonable approximation for Psoe{s,C() at sufficiently large values of s 
for all values of a 7^ 0. In this respect, the asymptotic behavior of the superstatistical NNS 
distribution is given by 

Psoe(s, a) - d exp (-C 2 s 2/3 ) , (29) 

where G\p are constants, unlike that of the NNS distribution obtained by Tsallis' non- 
extensive statistics [14], which asymptotically decays according to a power law. 

Figure 2 shows the evolution of PsoE(s,a) from a Wigner form towards a Poissonian 
shape as a decreases from 00 to 0. This distribution behaves similarly but not quite ex actly 
as any member of the large family of distributions. One of these is Brody's distribution [36], 
which is given by 

PiWs, 7) = a 7 * 7 exp (-a 7 ^ +1 /( 7 + 1)) , a 7 = -^—T^ +1 ( — !— ) . (30) 

7 + 1 \7 + 1/ 

This distribution is very popur but essentially lacks a theoretical foundation. It has been 
frequently used in the analysis of experiments and numerical experiments. The evolution of 
the Brody distribution during the stochastic transition is shown also in Fig. 2. The Brody 
distribution coincides with the Wigner distribution if 7 = 1 and with Poisson's if 7 = 0. 
On the other hand, the superstatistical distribution at a = is slightly different, especially 
near the origin. For example, one can use the small-argument expression of Mejer's G- 

function to show that lim Q > o, s >o -Psoe(s, a) = tt/2. In the mid-way of the stochastic 

transition, the agreement between the two distributions that is only qualitative. At small 
s, the superstatistical distribution increases linearly with s while the increase of the Brody 
distribution is faster. The large s behavior is different as follows from Eqs. (29) and (30). 
The difference between the two distributions decreases as they approach the terminal point 
in the transition to chaos where they both coincide with the Wigner distribution. 
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The superstatistical NNS distribution for systems in the midway of a stochastic transition 
weakly depends on the choice of the parameter distribution. To show this, we consider other 
;wo spa cing distributions, which have previously been obtained using other superstatistics 
lil [3 [ll The first is derived from the uniform distribution, considered in the original 



paper of Beck and Cohen |19j. The second is obtained for a ^-distribution of the parameter 
i], which is known to produce Tsallis' non-extensive theory. In the latter case, we qualify the 



NNS distribution by the parameter m = — j- — d — 2, where q is Tsallis' entropic index and 
d is the dimension of the Hamiltonian random matrix. This behavior is quite different from 
the conventional NNS which are frequently used in the analysis of experiments and nuclear 
experiments, namely Brody's and Izrailev's j^. The latter distribution is given by 

Pizraiicv(s, A) = As x exp (~^ 2 - j (B - A) sj, (31) 

where A and B are determined for the conditions of normalization and unit mean spacing. 
Figure 3 demonstrates the difference between the superstatistical and conventional distri- 
bution in the mid-way between the ordered and chaotic limits. The figures compares these 
distributions with parameters that produce equal second moments. The second moment of 
the Brody distribution is given by 

<* 2 >b d = 7 { ■ ( 32 ) 

We take 7 = 0.3, calculate (s 2 ) Brody and use the corresponding expressions for the second 
moment of the other distributions to find the value of their tuning parameters that makes 
them equal to (s 2 ) Brody . The comparison in Fig. 3 clearly shows that, while the considered 
three superstatistical distributions are quite similiar, they considerably differ from Brody's 
and Izrailev's distributions. 

The superstatistical distribution Psoe(s, a) can still be useful at least when Brody's 
distribution does not fit the data satisfactorily. As an example, we consider a numerical 
experiment by Gu et al. j3| on a random binary network. Impurity bonds are emp loyed to 
replace the bonds in an otherwise homogeneous network. The authors of Ref. 38} numeri- 
cally calculated more than 700 resonances for each sample. For each impurity concentration 
p, they considered 1000 samples with totally more than 700 000 levels computed. Their 
results for four values of concentration p are compared with both the Brody and supersta- 
tistical distribution in Figure 4. The high statistical significance of the data allows us to 
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assume the advantage of the superstatistical distribution for describing the results of this 
experiment. 



B. Unitary ensembles 

Now we calculate the superstatistical NNS distribution for a mixed system without time- 
reversal symmetry. Chaotic systems belonging this class are modeled by GUE for which the 
Wigner surmise reads 

32 f 4 \ 

P ^D) = —se W (-- -s*y (33) 

We again assume that the local mean spacing distribution f(D) is given by Eq. (13). The 
superstatistics generalization of this distribution is obtained by substituting (33) into (22), 



/ x 32a 3 2w™ ( s 4 « 2 2 

Psue(s, a) = nrV> ^ n/ — i - s 2 G™ a 3 + —^s 2 



1,-^0), (34) 



ir>L%G®(aP\0,i,l) Ud V ttD 2 

where D is given by (14). At large values of z, we use the large z asymptotic formula for 
the G-function to obtain 



32 6XP 
P SVE (s,a)^-s 2 — ^ 



3« ( a 3 /1 + ^ - 1 



(35) 



(1+421)" 

which clearly tends to the Wigner surmise for the GUE as a approaches infinity as in the 
case of a GOE. 

Figure 5 shows the behavior of Psue(s, ot) for different values of a ranging from to oo 
(the GUE). As in the case of the orthogonal universality, the superstatistical distribution 
is not exactly Poissonian when a = 0. Using the small argument behavior of Mejer's G- 
function, one obtains lim a >0s > P SOE (s,a) = 4/-7T. 



V. TWO-LEVEL CORRELATION FUNCTION 



The two-level correlation function is especially important for the statistical analysis of 
level spectra |4|. It is also directly related to other important statistical measures, such as 
the spectral rigidity A3 and level- number variance S 2 . These quantities characterize the 
long-range spectral correlations which have little influence on NNS distribution. 
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The two-level correlation function R 2 {E\, E 2 ) is obtained from the eigenvalue joint dis- 
tribution function P\ '(r), Ei, ■ ■ ■ , Em) by integrating over all eigenvalues except two. It is 
usually broken into a connected and disconnected parts. The disconnected part is a prod- 
uct of two level densities. On the unfolded spectra, the corresponding two-level correlation 
function can be written as Q, ^ \ 

X 2 (&, 6) = D 2 R 2 (D£ h (36) 

Here the disconnected part is simply unity and the connected one, known as the two-level 
cluster function, depends on the energy difference r = £1 — £2 because of the translation 
invariance. One thus writes 

X 2 (r) = 1-Y 2 (r). (37) 

The absence of all correlation in the spectra in the case of the Poisson regularity is formally 
expressed by setting all k-level cluster functions equal 0, and therefore 

X 2 Poisson (r) = 1. (38) 

We shall here consider the unitary class of symmetry. For a GUE, the two-level cluster 
function is given by 

U-W=(^)' (39) 
The two-level correlation function for mixed system described by the superstatistics formal- 
ism is given using Eqs. (7) and (26) by 

Xl VE (r) = jf ' f(D)±- 2 X^ { L )dD , (40) 

where we divide by (D~ 2 ) in order to get the correct asymptotic behavior of X 2 (r) — > 1 as 
r — > 00. Unfortunately, we were not able to evaluate this integral analytically in a closed 
form. The results of numerical calculation of Xf UE (r) for a = 0.5, 1 and 00 (the GUE) are 
given in Fig 6. 

VI. SUMMARY AND CONCLUSION 

We have constructed a superstatistical model that allows to describe systems with mixed 
regular-chaotic dynamics within the framework of RMT. The superstatistics arises out of 
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a superposition of two statistics, namely one described by the matrix-element distribution 
cxp [-r]Tr (H*H)] and another one by the probability distribution of the characteristic 
parameter 77. The latter defines the energy scale; it is proportional to the inverse square of the 
local mean spacing D of the eigenvalues. The proposed approach is different from the usual 
description of mixed systems, which model the dynamics by ensembles of deformed or banded 
random matrices. These approaches depend on the basis in which the matrix elements are 
evaluated. The superstatistical expressions depend on Tr(i^iJ) which is invariant under 
base transformation. The model represents the spectrum of a mixed system as consisting 
of an ensemble of sub-spectra to which are associated different values of the mean level 
spacing D. The departure of chaos is thus expressed by introducing correlations between 
the matrix elements of RMT. Spectral characteristics of mixed systems are obtained by 
integrating the respective quantities corresponding to chaotic systems over all values of D. 
In this way, one is able to obtain entirely new expressions for the NNS distributions and 
the two-level correlation functions for mixed systems. These expressions reduce to those of 
RMT in the absence of fluctuation of the parameter D, when the parameter distribution is 
reduced to a delta function. They can be used to reproduce experimental results for systems 
undergoing a transition from the statistics described by RMT towards the Poissonian level 
statistics, especially when conventional models fail. This has been illustrated by an analysis 
of a high-quality numerical experiment on the statistics of resonance spectra of disordered 
binary networks. 



VII. APPENDIX 



For sake of completeness, we give in this appendix the definition of the Meijer G-function 
as well as some of its properties, which have been used in the present paper. Meijer's 
G-function is defined by 





a>i, ■ 




(■ 






61, • 


••A ) 



=-f 

2niJ L 



nr =1 rfe + g )n; =1 r(i- aj - g ) 



z~ s ds, (41) 



where < n < p and < m < q while an empty product is interpreted as unity. The 
contour L is a loop beginning and ending at —00 and encircling all the poles of Y (bj + s) ,j — 
1, • • • , m once in the positive direction but none of the poles of Y (1 — a 3 — s) , j — 1, • • • , n. 
Various types of contours, existence conditions and properties of the G-function are given 
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Ill 



. The way by which integrals of the type considered in this paper are expressed in 



terms of the G-functions are described in 

The asymptotic behavior of Meijer's G-function, as \z\ — > oo, is given by |41| 
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(42) 



where a = q — p > 0, and a6 = |(1 — a) + J^| =1 &j — Sj=i a r m particular, the G-function 
that appears in this paper 



Gl>° 3 (z\b 1 ,b 2 ,h) = — 



1 



2ni J L T (1 - b x - s) T (1 - & 2 - a) V (1 - 6 3 - s) 
has the following asymptotic behavior 

Go'3 (z 6 2 , 6 3 ) ~ ^ z (bi+h+h-W exp (_ 3 *i/3) . 
v3 



z _s ds, (43) 



On the other hand, the small z behavior of Meijer's G-function 



42| is given by 



(44) 
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(45) 



Thus, the leading term in the expansion of G%° 3 (z 1 61, 6 2 , 63) in powers of z is given by 



g$ (z \h, b 2 , 63) « r (6 2 - &o r (63 - 61) A 

where 61 is the smallest of 6j. 

The implementation of Meijer's G-function in Mathematica 



(46) 



42 1 constitutes an additional 



utility for analytic manipulations and numerical computations involving this special function. 
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Figure Caption (Color on Line) 

Fig. 1. Level density for superstatistical orthogonal ensembles with parameter a = 0.2, 1 
and oo (the GOE limit). 

Fig. 2. Evolution of NNS distributions obtained by the superstatistics method for systems 
undergoing a transition from the GOE statistics to the Poissonian, compared with the 
Brody's distributions. 

Fog.3. Comparison between the superstatistical and conventional NNS distributions hav- 
ing equal second moments. 

Fig. 4. NNS distributions of geometrical resonances in random network, calculated by 
Gu et al. j3] compared with the Brody and superstatistical distributions. 

Fig. 5. NNS distributions obtained by the superstatistics method for systems undergoing 
a transition from the GUE statistics to the Poissonian 

Fig. 6. Two- level correlation functions obtained by the superstatistics method for systems 
undergoing a transition from the GUE statistics to the Poissonian. 
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